Phenomenological description of competing 
antiferromagnetism and d-wave superconductivity in high T c cuprates 



O 

o 

(N 

Of 
< 



CO 

c3 



c 

o 
o 



Bumsoo Kyung and A. M. -S. Tremblay 1 
Departement de physique and 
Centre de recherche sur les proprietes electroniques de materiaux avances. 
1 Institut canadien de recherches avancees 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1 
(April 23, 2002) 

In this paper the phase diagram of high T c cuprates is qualitatively studied in the context of com- 
peting orders: antiferromagnetism, d-wave superconductivity and d-density wave. Local correlation 
functions are estimated from a mean-field solution of the t — J Hamiltonian. With decreasing doping 
the superconducting mean-field and order parameter d begin to decrease below some character- 
istic doping x c ~ 0.2 where short-range antiferromagnetic correlations begin to develop. Dynamical 
properties that involve the energy spectrum, such as the normal state pseudogap, are calculated 
from effective interactions that are consistent with the above-mentioned local correlation functions. 
The total excitation gap A tg (in the superconducting state) and the normal state pseudogap A pg 
are in good agreement with experimental results. Properties of the condensate are estimated using 
an effective pairing interaction V e f / which takes into account (pair breaking) antiferromagnetic cor- 
relations. These condensate properties include condensation energy U(0), coherence gap A cff and 
critical field H C 2- The calculated coherence gap closely follows the doping dependence of T c or d, 



and is approximately given as A c 



A P9 within our numerical uncertainties. The system- 



atic decrease of superfluidity (d, {7(0), A cs , H C 2), and systematic increase of A pg and A tg with 
decreasing doping below x c have their natural explanation in our approach. The overall description 
is however qualitative since it does not appear possible to obtain results that are in quantitative 
agreement with experiment for all physical quantities. 

PACS numbers: 71.10.Fd, 71.27.+a 
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I. INTRODUCTION 

Fifteen years have passed since the discovery of high 
temperature superconductivity in Ln2- x Ba 2; Cu04 by 
Bednorz and Miiller, [jjj yet there is no widely accepted 
theoretical explanation of the many anomalous features 
found in these materials. This is due to a large extent 
to the lack of reliable theoretical methods to study this 
strongly-correlated electron problem. Dynamical Mean 
Field theory, for example, B works extremely well in 
three or higher dimensions but it lacks the momentum- 
dependent self-energy effects that are intrinsic to a two- 
dimensional situation with strong antiferromagnetic fluc- 
tuations. More recent dynamical cluster approaches cor- 
rect this defect but these methods are still in their infancy 
and are entirely numerical. || Q In this paper, we step 
back and take a more humble stand, seeking the simplest 
phenomenological approach that describes the main over- 
all features of high-temperature superconductors. Let us 
recall what these features are. 

First consider the generic phase diagram of a hole- 
doped cuprate L^-zSr^CuO^ in the doping (x = 1 — n) 
and temperature (T) plane. H| Near half-filling and at 
low temperature, antiferromagnetic (AF) long-range or- 
der appears and is destroyed by 2% doping concentration 
{x = 0.02). When x reaches 0.05, superconducting (SC) 



long-range order starts to appear and is also destroyed by 
about 30% doping. In between them, T c reaches a maxi- 
mum value at x ~ 0.15. The SC gap was found to have 
mainly d-wave character, in contrast to conventional 
BCS superconductors with an isotropic s-wave gap. 

Second, various recent experiments also show the ex- 
istence of a crossover temperature T* higher than T c in 
underdoped and optimally doped samples. Below this so- 
called pseudogap temperature T* , low- frequency spectral 
weight begins to be strongly suppressed in many phys- 
ical quantities. This anomalous phenomenon has been 
observed in angle resolved photoemission spectroscopy 
(ARPES), [0(l| specific heat, Jp) tunneling, @ NMR, 
|15| and optical conductivity. |16| Surprisingly the dop- 
ing dependences of T* and T c are completely different 
|l7],[L8| in spite of their possibly close relationship sug- 
gested by ARPES jll],(L|| and tunneling |L4| experiments. 
At optimal doping, where T c is maximum, various non- 
Fermi liquid (NFL) properties are observed in the normal 
state. These include the linear temperature dependence 
(up to 1000iv ) of ab-plane resistivity, the quadratic T 
dependence of Hall angle and so on. Far beyond optimal 
doping, there are indications that the normal state prop- 
erties are probably well described by the conventional 
Landau Fermi liquid. Near and below optimal doping, 
the superfluid density n s /m is systematically suppressed 
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with decreasing doping in spite of the increasing exci- 
tation gap amplitude, as measured by ARPES, in the 
superconducting state. The resulting 2A max / ksT c ratio 
is very different from the universal BCS value, 4.3 for d- 
wave symmetry. In the overdoped regime, however, the 
SC properties are probably well explained in the conven- 
tional weak-coupling BCS theory, although this has not 
yet been checked in detail. 19 

The approach that we take is as follows. First, we 
do a simple mean-field analysis of the t — J Hamilto- 
nian whose results are used only to obtain estimates 
for zero-temperature order parameters and to compute 
equal-time local correlation functions (such as double oc- 
cupancy). Dynamical quantities, such as the response 
functions and single-particle spectral weight, are com- 
puted by using the values of the local correlation func- 
tions obtained in mean field to fix the value of the renor- 
malized interaction vertices. We do not consider the pos- 
sibility of spin-charge separation |2^,^l[ but we take into 
account the no-double-occupancy constraint by replacing 
projectors in the t — J Hamiltonian by their mean-field 
bosonic value, namely doping x. This is an effective-mass 
approximation. We cannot address the issue of a quan- 
tum critical point (QCP) 0] at x c even if it exists at 
the mean-field level. Non-Fermi liquid behavior above 
the pseudogap temperature is also beyond the domain of 
validity of our approach. The regime we consider is one 
where the temperature would normally be low enough for 
a Fermi liquid to appear but where fluctuations destroy 
that phase, especially in the underdoped regime. 

The paper is organized as follows: In Section [n] the 
t — J Hamiltonian is decoupled, in a mean-field approx- 
imation, with competing (SC, AF and d— density wave) 
order parameters. The mean-field Hamiltonian is diag- 
onalized to compute mean-field ordering temperatures 
T^ IF , T^ IF and mean-field order parameters d and m. 
Dynamical quantities such as the normal-state pseudogap 
are computed and compared with experimental results in 
Section [II . In Section |^ the results of the previous sec- 
tions are gathered to discuss the phase diagram of high T c 
cuprates. In Section ^ we discuss some of the properties 
of the SC condensate, such as the condensation energy, 
upper critical magnetic field and the coherence gap, using 
a model that takes into account the effect of (pair break- 
ing) AF correlations in an effective manner. The present 
results are compared with some of the leading theories 
for the high T c superconductivity in Section ^ as well 
as with experiments. The last Section summarizes the 
results and discusses its main limitation. 



II. MEAN-FIELD T - J HAMILTONIAN 

Right after the discovery of high temperature super- 
conductors, Anderson [^3| first proposed the one-band 
Hubbard model as the simplest Hamiltonian that might 



capture the correct low energy physics of copper oxides. 
The t — J model is known to be the large U limit of the 
Hubbard Hamiltonian under certain assumptions. The 
t — J model is described by the Hamiltonian 



(1) 



where Ci i(T destroys an electron at site i with spin a on 
a two-dimensional square lattice, t is the hopping ma- 
trix element between nearest neighbors < i,j > and J 
denotes superexchange coupling. Double occupancy of 
two electrons at the same lattice site is forbidden in the 
Hilbert space corresponding to Eq. (Jl|) as can be seen by 
the projection operator (1 — n^_ CT ) in the hopping term. 
The chemical potential [i controls the electron density n. 
Finally, Si and rij are spin and charge density operators, 
respectively, and they are defined by 



a,/3 



C i,cr Ci ' a ' 



(2) 



where the components of a are 2x2 spin Pauli matrices. 

In general the (static) mean-field description is not ac- 
curate for strongly correlated electron systems such as 
the Hubbard and t — J models in the physically rele- 
vant regime. Nevertheless it is known that, for example, 
the mean-field description of the antiferromagnetic state 
at x = is quite accurate, even at large U. For exam- 
ple, spin waves with accurate spin-wave velocities can be 
obtained by performing RPA in the ordered phase. |24|] 
We take the point of view that mean-field theory can 
be useful to select the leading correlations when there 
are several competing order parameters in the low energy 
sector. Due to the neglect of spatial and quantum fluctu- 
ations in mean-field theory, some caution is necessary in 
interpreting the results. First, the mean-field phase line 
is to be decreased in temperature to take into account 
the effect of quantum fluctuations (such as Kanamori- 
Brueckner screening p5[). Second, since we work in two 
dimensions where antiferromagnetic long-range order is 
prohibited by the Mermin- Wagner theorem, the renor- 
malized mean-field Neel transition line has to be inter- 
preted as the onset of the corresponding short-range cor- 
relations instead of as a true thermodynamic phase tran- 
sition. Whether given short-range correlations eventu- 
ally grow to long-range order at low temperature or not, 
should be answered by studying how the zero-frequency 
correlation length grows when temperature is lowered. 
This question, as well as that of dynamical or quantum 
fluctuations that may take place on very short distance 



2 



scales S can be examined by the fluctuation theory that 



we describe in Section [II 



In the spirit of a mean-field approximation, terms with 
more than two operators should be decoupled in all pos- 
sible ways. In principle there are infinitely many ways of 
decoupling the Heisenberg part of the Hamiltonian. In 
this situation, guidance from experiments and past the- 
oretical considerations is helpful to find the most impor- 
tant leading correlations of the model, correlations whose 
stability is studied a posteriori. We consider three com- 
peting mean-field order parameters, m, d, and y for AF, 
d-wave SC, and d-density wave [^6j orders respectively, 
which are defined |27|] by 



n=(m) = l/(2iV)$><4 +4ff%r ), (3) 

k.a 

d = (a) = l/NY,<t>d(k){% r c_ u ) , (4) 

y = (y) = if (2JV) ]T MW M<r %J, (5) 



k.a 



where N is the total number of lattice sites, (j) d (k) — 
cos k x — cos k y is the d — wave form factor and Q is 
the (commensurate) AF wave vector (71", it) in two dimen- 
sions. The extended s— wave analogs of d— wave super- 
conductivity and d— density wave orders are also present 
in the t—J Hamiltonian, but they are always less relevant 
in mean-field studies. In this paper we restrict ourselves 
to a uniform solution. At the end of the paper the issue 
of inhomogeneous modulation of spin and charge degrees 
of freedom will be briefly discussed in the context of the 
present results. A similar mean-field decoupling for the 
AF and SC channels only was previously considered by 
several groups. |28|j29| l In the present study the projected 
operator Cj i<T (l— nj- a ) is separated into the original elec- 
tron operator Cj )(T times the expectation of the bosonic 
operator (1 — nj t -„) which is taken as 2 1 / 2 . As long as 
the no-double-occupancy constraint is globally imposed, 
the average number of electron n e in our study is equal to 
1 — nh with x — rih- Clearly, at half-filling, one recovers 
the Heisenberg Hamiltonian. In addition, as we shall see 
later, for T = 0, n = 1, the no-double-occupancy con- 
straint is satisfied exactly because the AF solution allows 
only one electron per site with the full moment. This 
would not be the case if we considered only d— wave and 
d— density wave mean-field solutions. 

For a more systematic implementation of this con- 
straint beyond the mean-field level, the slave boson rep- 
resentation can be useful. |5(|[3l| In the slave boson rep- 
resentation, an electron is explicitly decomposed into a 
spinon (fermion) and a holon (boson), a — f^ a bi. The 
latter bosonic operator corresponds to (1 — ni_ CT ) in our 
approach and our global constraint 1 = n e + corre- 
sponds to the no double-occupancy constraint in slave 
boson mean-field theories. 



In terms of mean-field order parameters m, d, and y, 
the mean-field t—J Hamiltonian reads 

Hmf = e(k)ct ct — AJmih— Jdid + $) — Jyy + F Q , 

(6) 



k . o 



where 



.Fq = N(2Jm 2 + Jd 2 + ^y 2 - fi) 



(7) 



e(k) ~ — 2te(cos 

kx ~\~ cos ky) /1 with x the hole density. 
The correlated hopping is taken into account through x 
in a mean-field spirit, as explained above. In other words, 
the no-double-occupancy constraint manifests itself here 
as an effective mass, not as a change in the number of 
carriers. One has to keep in mind however that, at this 
level, e(k) and the corresponding effective mass are only 
bare values that are renormalized by interactions. This 
is discussed further at the end of the following section. 

Previous studies have suggested that the d— density 
wave order parameter becomes important only quite close 
to half-filling. |}2| Here we find that this order parameter 
is never a self-consistent solution of the mean-field equa- 
tions when mean-field AF order is explicitly considered 
in Eq.(^|). Hence, in the following, we restrict ourselves 
to AF and SC order parameters. By introducing a four 
component field operator 

k 

Eq. (|^) may be written in a more compact form to reflect 
the competition between AF and SC order, 



H 



AIF 



F 



The prime symbol on the summation restricts the sum- 
mation over wave vectors to the magnetic Brillouin zone 
in order to take into account the doubling of the unit cell 
in the presence of (commensurate) AF order. Small in- 
commensuration would not change appreciably the value 
of the local correlation functions that we will compute 
with the mean-field solution, hence we limit ourselves to 
the commensurate case. The matrix in the last equa- 
tion is given by 



M k 



( e{k) -Jd(f) d {k) -2Jm 

-Jd<j) d {k) -e(k) 

-2Jm e{k + Q) 

\ -2Jto Jd(j) d {k) 



\ 





-2Jm 
Jd(f> d (k) 
-e{k + Q) j 



(8) 



The energy eigenvalues of yield four branches to the 
energy dispersion ±E±(k) 
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±90) ] 1/2 , 



{Jd<t> d {k)f 



(9) 



where g(k) is given as 

g(k) = [{e\ - 4 +Q -) 2 /4 + ((eg + e £+(3 )(2 Jm)] 1/2 . (10) 

The free energy is easily obtained either from the trace 
formula or from the Feynman theorem 



F 



-2T 



k a=± 



F 



Two mean-field equations are obtained from the station- 
ary condition on F with respect to the corresponding 
order parameters, = ^ = 0, and one more unknown 
constant fi is determined by the thermodynamic relation 
n = 1 — x = —^f-- The resulting three equations arc 



m) + a- 



t a=± 



2g{k) 



1 . ,,§EjJt) 
x — tanh( — ) , 



E a {k) 



d = ^EE <p 2 d (k)(Jd)-^ tanh( 



5 a=± 



+a- 



+a- 



2g(k) 

2{2Jmf{ej: + e- k+Q ), i 



E a {k) 



tanh( 



pE a {k) 



)• 



(11) 



The spin-triplet order parameter (cj + q |C_£ j_) is "dy- 
namically generated" |?3| when the two mean-field order 
parameters to and d coexist, even when it is not explic- 
itly included in the mean-field factorization. We have 
checked that including this order parameter in the mean- 
field factorization changes the value of local correlation 
functions, such as (|A,j(0)| 2 ) which we need below for the 
dynamical calculation, by less than one percent. The val- 
ues of to and d, on the other hand, change by about 10%. 
The inclusion of this order parameter would then change 
the results of the present paper only quantitatively and 
at the few percent level only. We do not include it for 
now for simplicity. It should be included at a later stage 
when refining the approach. 

Even though it is only an intermediate step in our 
calculations, for completeness we show in Fig.[l] the cal- 
culated mean-field phase diagram and zero-temperature 




FIG. 1. (a) Calculated mean- field phase diagram in doping 
(x = 1 — n) and temperature (T) plane and (b) mean-field or- 
der parameters for t/J = 3. T^ F and T C MF are mean-field AF 
and SC ordering temperatures, while m and d are mean-field 
AF and SC order parameters. T^f F (m = 0) and d(m = 0) 
(dashed curves) are SC ordering temperature and order pa- 
rameter, respectively, when the mean-field factorization of the 
interaction term (J) is in the SC channel only. The dashed 
curve for Tff F (d = 0) is also shown in (b) but it is indistin- 
guishible from the case where d takes its mean-field value. 



mean- field order parameters to and d for t/J = 3. The 
AF order parameter dominates near half-filling and van- 
ishes beyond x — x c ~ 0.20. On the other hand, d- 
wave SC correlations keep growing with decreasing dop- 
ing down to x c ~ 0.20 where the AF order parameter is 
non-zero. Below x c , we observe that T C MF and d begin 
to monotonically decrease with decreasing doping. The 
physical origin of this decrease of T^ F and d should be a 
combination of the effects of local moment formation and 
of the onset of AF correlations that break time-reversal 
symmetry, both of which lead to the destruction of the 
SC state. The importance of AF correlations in reduc- 
ing T c to the experimental value (~ 100K) was recently 
pointed out by Weng. §34| In between half-filling and x c , 
AF and SC orders coexist at the mean-field level. The 
overall features are similar to previous results Q based 
on a phenomenological model in which SC and AF orders 
come from different interaction terms. 

The dashed curves in Fig.|l|(a) and (b) are T^ IF and 
d when the interaction piece of the Hamiltonian is de- 
coupled only in the (spinon) pairing channel in the slave 
boson mean-field theory ]35[ | or, equivalently, when m is 
forced to zero in the present approach. This shows that, 
by itself, the diverging effective mass is not enough to re- 
duce T^ F and d to zero at half-filling. In the slave-boson 
theory, T c MF and d keep increasing with decreasing dop- 
ing. T^ IF is interpreted as pseudogap temperature T* 
and the corresponding ground state is dubbed as the res- 
onating valence bond (RVB) state. ]23| In this theory, the 
true T c is governed by the Bose-condensation tempera- 
ture of holons, which scales as x. Then the maximum 
T c occurs at very small doping concentration (x < 0.05). 
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More recent versions of this slave-boson mean-field the- 
ory seem to correct this defect [ [36| . 

In the present study, the behavior of the mean-field or- 
der parameters reflects the competition between AF and 
SC correlations. As can be seen in Fig.[l], SC correlations 
are most significantly modified due to the presence of 
(pair breaking) AF correlations, while the latter are vir- 
tually unchanged by the former. When Eq. ||is converted 
into an effective Hamiltonian where the competing order 
parameters are expressed in terms of operators we find 



H = 



2Jrhrh - J$d- i Jyy. (12) 



In this notation, it is clear that interactions that gen- 
erate AF are stronger than those that generate SC by 
a factor of two and that the interactions that generate 
d— density wave order are even smaller. Note that to be 
able to make this comparison, the order parameters are 
normalized j2?J such that the largest value that either 
one of them can take in a mean- field solution is near 1/2 
(when all other order parameters are forced to zero). In 
particular then, even if we drop the d— density wave, this 
mean-field t— J Hamiltonian is not SO(5) symmetric ]37j 
since in that case, the mm and d^d terms would have the 
same coefficients, so as to form a five-dimensional super- 
spin vector. [B8J 



III. DYNAMICAL PROPERTIES 
PSEUDOGAPS 



EXCITATION 



In a previous study of the two-dimensional attractive 
Hubbard model |39| it was shown, by comparisons with 
Quantum Monte Carlo calculations, that the BCS mean- 
field ground state gave a good estimate (few percent) of 
double-occupancy in the normal state when the interac- 
tion strength satisfied \U\ > 3t and when the temper- 
ature is reasonably low (T of order 0.5i or less). We 
expect this should remain true for correlation functions 
that are local in space and time. On the other hand, for 
dynamical quantities mean-field theory itself is clearly a 
bad estimate. Instead of a true antiferromagnetic gap as 
in Fig.0, we expect to have a pseudogap. 

In Ref. (40|, the procedure used to compute the pseu- 
do gap size A^g was as follows. Let A^(z) be defined by 



where g(6) is a d-wave form factor given, in real space 
representation, by g(S) = 1/2 for 5 = (±1,0), —1/2 
for 5 = (0,±1), and otherwise. The <i-wave symme- 
try spin-singlet local correlation function, ( | A.^ (0) | 2 ) , is 
computed from the ground state mean-field solution ]ic|] 
but with only m different from zero. In other words, it 
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FIG. 2. (a) Calculated pseudogap size at T — 0, A ps , 
denoted as filled diamonds for t/J = 3.0 and J = 125 
meV. The open circles, diamonds, squares, up-triangles, and 
down-triangles are the pseudogap size extracted from suscep- 
tibility, heat capacity, ARPES, NMR, and resistivity measure- 
ments, respectively, by Tallon and Loram. [^| The dashed 
curve denotes (spinon) pseudogap size predicted by the slave 
boson mean-field theory of the t — J model, (b) Calculated 
total excitation gap at T = 0, At g (empty diamonds). The 
empty circles with error bars denote ARPES leading edge gap 
by Ding et al. Q 



contains only AF correlations. This correlation function, 
(|A<j(0)| 2 ), is then used in the fluctuation-dissipation the- 
orem to obtain a renormalized vertex that allows one to 
estimate dynamical susceptibilities in any channel. 
The latter come in the computation of both the two- 
particle and one-particle quantities. An example of such 
a calculation is given in section^. The pseudogap ap- 
pears in the renormalized classical regime j42| of the sin- 
glet d-w&ve fluctuations in a region of the Fermi surface, 
namely the zone edges, where the thermal de Broglie 
wavelength is so small that it can become smaller than 
the ci-wave zero-frequency correlation length even if the 
latter is of order only one lattice spacing. To be more 
specific, for x = 0.1 for example, the pseudogap appears 
near the zone edge when T = 0.40 J. At this temperature, 
the characteristic energy of the d— wave singlet fluctua- 
tions is 0.31J while the corresponding correlation length 
estimated from the wave vector dependence of the u> = 
d— wave singlet susceptibility around q = is 0.48 lattice 
spacing. 

In Fig. |](a) the calculated pseudogap size A pg (filled 
diamonds) from Ref. p(| is plotted along with the pseu- 
dogap energy (empty symbols) extracted from various ex- 
periments by Tallon and Loram. pl| Note that the pseu- 
dogap rapidly takes its asymptotic low-temperature value 
when it opens up at finite temperature. That is the rea- 
son why we plot the value it would have at T = if the su- 
perconducting phase did not appear. Since in the above 
calculational procedure, the local d— wave singlet correla- 
tion function (|Ad(0)| 2 ) that is used to compute the dy- 
namical singlet susceptibility is computed from static an- 
tiferromagnetic mean-field order, the corresponding nor- 
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mal state pseudogap A pg vanishes at x ~ x c . The cal- 
culated A pg is in good agreement with experimental re- 
sults. [Q Through interlayer tunneling spectroscopy of 
Bi2Sr2CaCu20s+s Krasnov also observed that the 
pseudogap decreases approximately linearly with doping 
and vanishes at the critical doping, x c ~ 0.19. 

The dashed curve is the pseudogap size predicted from 
the slave-boson mean-field theory of the t — J Hamilto- 
nian. [||| In fact the (spinon) pseudogap size obtained 
in the latter theory is the same as 2Jd calculated from 
Eq.(^|) with m — 0. It is much larger than the experimen- 
tally extracted pseudogap energy throughout the inter- 
esting doping range. Furthermore there is no indication 
that it vanishes near x c , in conflict with experimental 
results. Although spin degrees of freedom are the main 
driving force for the pseudogap behavior in the above 
theory and in ours, the detailed implementation is quali- 
tatively different. In the slave-boson mean- field theory a 
pseudogap is caused by a pairing among spinous, while 
in our theory it is short-range spin correlations of elec- 
trons that induce an effective attraction in the particle- 
particle channel. The spinon pseudogap in the former 
theory has nothing to do with local moment formation 
and their short-range AF correlations with decreasing 
temperature. 

The total excitation gap A tg , at T = 0, is obtained 
by allowing both d and m to be different from zero in the 
calculation of the quantity (|Ad(0)| 2 ) that is used in the 
fluctuation-dissipation theorem to obtain the renormal- 
ized vertices entering the calculation of dynamical quan- 
tities. The calculated total excitation gap (or ARPES 
leading edge gap or SC gap), A tg , at T — is also plot- 
ted in Fig. ||(b) as a function of doping with t/J = 3. 
The open circles with error bars are experimentally de- 
termined leading edge gap by Ding et al. E| in their 
ARPES measurement. Except a few points, which may 
be due to the effect of impurities as noted by the authors, 
our calculated A tg is in reasonable agreement with exper- 
imental results and decreases more or less linearly with 
doping. In addition, the velocity defined by 

A tg . k F 
v 2 = —7= sin — = 
s/2 V2 

ranges from 1.0 to 2.0 x 10 6 cm/s when x varies from 0.2 
to 0.1. Experiments on various compounds and dopings 
in that range give results of order 1.0 x 10 6 cm/ s accord- 
ing to Ref. @ More recent ARPES results Bi-2212 @ 
give a ratio of V2I A tg that can decrease by about a factor 
two from optimal to underdoping, while transport mea- 
surements suggest a much stronger doping dependence. 
p9[ Note that in our calculations of A pg and A tg we 
use the experimentally determined value of exchange cou- 
pling J — 125 meV. The only adjustable parameter is the 
ratio t/J that we take, as in the previous sections, equal 
to 3. 



A tg is always larger than A pg due to the additional 
contribution from d 7^ jl0| to the local spin-singlet am- 
plitude (|Arf(0)| 2 ). Since the SC order parameter van- 
ishes at T c , the SC gap below T c continuously evolves 
into the normal state pseudogap above T c with the same 
momentum dependence and magnitude. 

Finally we comment on the Fermi velocity in the nodal 
direction. |l7j Tabulation of angle-resolved photoemis- 
sion (ARPES) data for the Fermi velocity [^6| along the 
zone diagonal in the range 0.1 < x < 0.15 lead to Fermi 
velocities in the range range from 1 to 2.5 x I0 7 cm/s. 
However, recent ARPES data on LSCO seem to sug- 
gest that the Fermi velocity along the diagonal is dop- 
ing independent and of order 2.5 x 10 7 cm/s. A weak 
doping dependence for the limiting low frequency Fermi 
velocity is seen in Bi-2212. |0 Our bare Fermi velocity 
varies linearly with x, but to compare with ARPES ex- 
periment, one should include the effect of residual inter- 
actions. The superconducting d-wave part of the inter- 
action does not renormalize the Fermi velocity along the 
nodal direction, but antifcrromagnetic fluctuations do. 
For x = 0.15 and T — 0.05 J, we find that the bare veloc- 
ity, estimated from the dispersion of the spectral weight, 
renormalizes from 1.27aJ to 1.04 aJ. For x = 0.1, the 
peak in the spectral function is very broad and in fact 
the quasiparticle picture does not strictly applies since 
dYj/doj > 0. Nevertheless, if we do like in experiment 
and measure the dispersion of the spectral weight maxi- 
mum, we find that for T between 0.2 J and 0.1J, the 
Fermi velocity renormalizes up (since dYi/duj > 0)from 
0.85 aJ to 0.95 ± 0.05 aJ. Hence, the physical value 
of the velocity is much more doping independent than 
the bare value suggests. Such doping-independent val- 
ues were also found in variational calculations that take 
into account the no-double occupancy constraint. l5l| For 
t/J = 3, and J = 125meV, our Fermi velocity in physical 
units is thus about 0.7 x 10 7 cm/s, which is smaller than 
experimental values by a factor 3 to 4. 

IV. QUALITATIVE PHASE DIAGRAM OF HIGH 
T c CUPRATES 

Applying the fluctuation approach of the previous sec- 
tion with the effective theory of the next section, one can 
use the resulting dynamical SC susceptibility to estimate 
the value of the correlation length £ sc at the tempera- 
ture where order appears according to mean-field the- 
ory. This length, £ sc , is of order one lattice spacing. We 
know that quantum renormalization effects (Kanamori- 
Brueckner type) should decrease even the mean-field T c 
to a lower value T C MF . In addition, thermal fluctuations 
decrease it even further. If we look empirically at the 
value of ^ sc when one reaches the experimentally deter- 
mined superconducting transition temperature T c , one 
finds that at that temperature, £ sc ~ 3 — 4. We know 
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FIG. 3. Calculated phase diagram from fluctuation theory. 
T c (empty circles) is the SC ordering temperature determined 
by fluctuation theory in which the correlation length £ sc be- 
comes 4 lattice spacing. The solid curve is an interpolation of 
these points given as T c = -6122(x - 0.20) 2 + 138 (K). T C MF 
(dashed curve) and T* (filled diamonds) denote mean-field 
SC ordering temperature and pseudogap temperature, respec- 
tively. 



that when the in-plane £ sc is sufficiently large, it should 
induce a three-dimensional SC transition. Since £ sc ex- 
ponentially increases at low temperature in the renor- 
malized classical regime, [^5[^9) the temperature where 
£ sc ~ 3 — 4 is not much different from that where it is 
much larger. The empty circles in Fig. |^ are the value 
of T c estimated from the above procedure (£ sc = 4). 
Unlike T^ IF , T c looks like a parabola centered at x c . 
The solid curve is an interpolation formula for the data 
points given as T c = -4.22(a; - 0.20) 2 + 0.095 (units of 
J) or T c = -6122(a; - 0.20) 2 + 138 (K). In the same plot, 
the dashed curve is the mean-field SC temperature T C MF 
and the filled diamonds are the pseudogap temperature 
T*. One expects that 0(2) SC fluctuations associated 
with superfluid stiffness should come into play only for 
t 'mf >T>T C , as found by Corson et al. §l| in their 
terahertz spectroscopy measurements. This result con- 
trasts with the SC fluctuation scenario for the pseudogap 
where SC fluctuations appear for T* > T > T c . The crit- 
ical region for SC fluctuations appears larger at under- 
doping than at ovcrdoping because of its small superfluid 
density. A rough estimate for the superfluid density is 
n/m* which scales like x (1 — x) in our approach because 
n *~ (1 — x) and 1/m* ~ x. The pseudogap temperature 
T* falls from a high value onto the T c line |4(J instead 
of sharing a common line with T c in the overdoped re- 
gion. This is because in our approach the d-wave singlet 
fluctuations that lead to the normal state pseudogap T* 
are induced by short-range spin correlations. The latter 
compete in the low energy sector with SC correlations to 
make T c go to zero near half-filling. |53[] Note however 



that since O (2) SC fluctuations also cause a pseudogap 
close to T c , (3^] some experiments may suggest that the 
T* line continues smoothly onto the dashed line (or more 
precisely the T C MF line) in the overdoped regime. 

In the present study T c and T^ F have their maximum 
at x c ~ 0.20 where condensation energy U(0), mean-field 
order parameter d and the coherence gap A cff discussed 
in the following section are also maximized. In high T c 
cuprates, the highest T c is at x = 0.16 slightly lower than 
the critical doping [x = 0.19) where the strength of su- 
perconductivity is maximum. j fl8[ This subtle difference 
is beyond the scope of the present study. 

Although short-range AF correlations persist up to x c , 
they don't lead to long-range (commensurate) AF order, 
remaining only as short-range order even at low temper- 
ature, unless the electron density is close to half-filling. 
For instance, at x = 0.15, the AF correlation length £ y 
computed from fluctuation theory saturates at ~ 4 
at the lowest temperature studied while the SC correla- 
tion length £ sc exponentially diverges. Thus the long- 
range (commensurate) AF phase boundary is not shown 
in Fig. |[ It is expected to lie close to half-filling where 
the present formulation is not valid for AF spin fluctua- 
tions since the effective bandwidth, 8xt becomes smaller 
than interaction strength. EG] 



V. THE CONDENSATE 

At finite temperature algebraic (Kosterlitz-Thouless) 
superconducting order can develop in two dimensions, 
contrary to AF order, which is prohibited by the Mermin- 
Wagner theorem. Since short-range AF correlations 
can produce a (dynamical) pseudogap on short-distance 
scales while thermodynamic quantities, such as the SC 
condensate, are calculated in the static and long-distance 
limit, wc should integrate out AF correlations to study 
the low energy physics of the SC state while taking into 
account the presence of the AF correlations in an effec- 
tive manner. In principle, one should be able to do this 
and predict the onset of superconductivity from the in- 
teracting Green function of the previous sections, which 
exhibits a pseudogap. However, to do this calculation, 
one needs the corresponding irreducible vertex. It is not 
strictly correct to use the Thouless criterion in the form 
1 — JGG = to find the superconducting T c because one 
cannot use the bare vertex J to do a calculation with 
dressed Green functions. Since it is not known yet |25| ] 
how to obtain a reliable approximation for irreducible 
vertices in the pseudogap regime, we have to proceed 
otherwise. Calculations with constant renormalization to 
the vertices and bare Green functions are expected to be 
more reliable. The phenomenological correlation length 
criterion of the previous section is such a calculation. 

We want to take into account the presence of short- 
range AF correlations. To proceed with a calculation 
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(a) t/J=3 




1-n x = 1- 

FIG. 4. (a) Effective strength of pairing interactions in the 
presence of (pair breaking) AF correlations and (b) calculated 
condensation energy at T = 0, as a function of doping con- 
centration. 

that is done with constant renormalization to the ver- 
tices and bare Green functions, we have to make an ad- 
ditional hypothesis. We construct an effective low-energy 
Hamiltonian 



H, 



- V t 



3// ^A+(i)A d (i) (13) 



such that the effective attraction V e ff leads, at T = 0, 
to the same value of the SC order parameter as that ob- 
tained from the full mean- field theory (Fig.^b)). We 
will use this effective Hamiltonian only to compute prop- 
erties of the condensate, including £ sc of the previous 
section. As long as superconductivity, or condensation, 
is concerned, the above Hamiltonian already takes into 
account the effect of (pair breaking) short-range AF cor- 
relations in an effective manner. 

Above x c , V e f f is given by J, because m vanishes there. 
We take the vanishing of m as signaling that the antifcr- 
romagnetic correlation length has become smaller than 
one lattice spacing. Below x c , V e ff rapidly vanishes with 
decreasing doping in Fig.[|(a). As a first application of 
the above effective Hamiltonian, let us calculate the con- 
densation energy at T = 0. The total energies in the 
SC state and in the normal state are, in a mean-field 
approximation, given by 

E d = (y d \H eff -fiN\y d ) 

r. k 



E n = {V n \H eff - (iN\y n ) 

= EK £ *-'*)-lte-M)|], 



(14) 



where — 



(e- k -^ + (dMk)V e ff) 2 



The T = 

condensation energy £7(0) = —(Ed — E n ) is plotted in 
Fig. ^b) as a function of doping concentration with 
J = 125 meV. J/(0) is maximum at x c , and rapidly 



decreases beyond and below x c . This is because U(0) 
scales like d 2 . This feature is qualitatively consistent with 
specific-heat measurements by Loram et al. |E3 The max- 
imum U(0) in our calculation is about 14A'/unit cell at 
x c , two to three times larger than what the above authors 
obtained (5 — 6AT/unit cell). It is clear that the actual 
condensation energy should be smaller than what we cal- 
culate since fluctuation effects are completely neglected 
in the above calculation. If we had stuck to a mean- 
field calculation with both AF and SC order present, we 
would have estimated the condensation energy from the 
difference between the ground-state energies computed 
with d finite and with d forced to zero. For x > x c we 
would have found the same result as on Fig. ^(b) and for 
x < x c we would have also found a result that decreases 
with doping, but it would have been larger than what is 
found above, suggesting that the above effective Hamilto- 
nian approach is indeed a better way to take into account 
T = long-range SC order in the presence of short-range 
AF correlations. 

Another quantity of interest is the coherence-gap A cg . 
Deutscher |35| has proposed that this gap, accessi- 
ble through Andreev reflection experiments with point- 
contact spectroscopy, can be different from the single- 
particle gap observed in tunneling. Recent experiments 
p6| seem to confirm this. To estimate this coherence gap, 
one can use the mean-field solution of the effective BCS 
model Eq.()l3|). However, a better estimate can be ob- 
tained by computing the precursor pseudogap with the 
method explained below. Indeed, it is known ]3S| ] in the 
attractive Hubbard model that this precursor pseudo- 
gap reaches its asymptotic low-temperature value very 
rapidly. In addition, this precursor pseudogap contains 
quantum renormalization effects that are absent from a 
pure mean-field calculation. The following calculation is 
in the spirit of Ref. First we compute how much 

the d-wave pair correlation function Xpp is enhanced over 
Xp P by applying the local-pair sum rule for \vv 



q 



= <|A d (0)| 2 ) . 



(15) 



The right-hand side of the above equation is evaluated in 
the mean-field state of Eq. |l^ and \pp is related to \p P 
through the rcnormalized vertex V pp 



x PP (q) 



x pp(q) 



1 - V pp x%(q) ' 
where the irreducible susceptibility is defined as 



x° P p(q) 



-E( 



t>d(k) + Mq- k)) 2 G°(g - k)G°(k) . 



Finally the following self-energy is used to estimate the 
effect of d-wave pairing correlations on quasiparticles. 
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FIG. 5. (a) Calculated total excitation gap At g (empty di- 
amonds), pseudogap A pg (filled diamonds) and condensation 
gap A cg (empty circles) at T — as a function of doping. 
The dashed and long-dashed curves denote Atg — A pg and 
Ap g , respectively, (b) The total excitation gap A 



A? 

tg 



I 9 



as a function of temperature for three representative doping 
levels (x = 0.08 for underdoping, x = 0.15 for optimal doping 
and x = 0.22 for overdoping). 



1 rp 

Spp(fc) = -jV e ffV pp — 2J 



(Mk) + Mq~ k)) 2 Xpp (q)G°(q - k) . (16) 

This procedure is analogous to that used earlier in the 
case of the attractive Hubbard model (s-wave symmetry) 
p9[ to obtain good agreement with QMC calculations. 
Even though the above formula has not been obtained 
with the same rigor, one expects it to give a good esti- 
mate of what happens in the case of <i-wave symmetry. 

Figure |5|(a) shows A cg (empty circles) together with 
A tg (empty diamonds) and A pg (filled diamonds) as a 
function of doping at T = 0. Below x c , A cg monoton- 
ically increases with doping and above x c it reduces to 
At s . The doping dependence of A cg is qualitatively sim- 
ilar to T C MF and d, as found in the experiments of Ref. 
|p6| . It is surprising to observe that A cg , A tg and A pg 
are more or less linear in doping concentration, in spite 
of the fact that they are obtained through highly non- 
linear equations. In fact the linear behavior of A tg and 



When A tg ~ A pg + A cg is generalized to finite temper- 
ature, the total excitation gap is composed of the nor- 
mal state pseudogap and the coherence-gap in Fig.||(b). 
A ps is virtually temperature independent from T* to 
T = 0, because it is caused by short-range spin corre- 
lations of order one lattice spacing. On the other hand, 
we expect that A cg is strongly temperature dependent 
since it opens up at T c . For purposes of illustration, we 
added a BCS-like temperature dependence to A cg . Its 
zero temperature- value and T c were estimated above. As 
one can see from |5|(b), at large underdoping (x — 0.08) 
A tg (or ARPES leading edge gap) is mainly given by 



= 0.15) it is more or 
below T c . At large 
thus leading to the 
conventional BCS behavior. This feature is also qualita- 
tively consistent with ARPES experiments. flTl]jL2f 



Apg, while at optimal doping [x 
less equally made of A pg and A c 
overdoping (x = 0.22) A 



tg 



A pg in doping are observed by Tallon and Loram [[18] 

and by Ding el al. Eq] in their respective experiments. the strong temperature dependence of the coherence gap 



A similar behavior of A tg is used by Chen et al. |57| to 
calculate the temperature dependence of some physical 
quantities within their SC fluctuation theory of the pseu- 
dogap. The crucial difference between their theory and 
ours is that their A pg vanishes at T = so that when 
superconductivity is completely destroyed, for instance, 
by a strong magnetic field, their ground state becomes a 
perfect Fermi liquid. In our picture, however, we predict 
that the ground state has a pseudogap. Q Our calcu- 
lated temperature dependence of the total excitation gap 
(or ARPES leading edge gap from the experimental point 
of view) contrasts with some theories based on the idea 
of spin-charge separation. In these theories the only dif- 
ference between A(k,w) above and below T c is the line 
shape due to the Bose condensation of holons, namely, 
the appearance of the quasiparticle peak below T c with- 
out change of the gap magnitude itself. At large under- 
doping, there is indeed not much change in the ARPES 
leading edge gap [ ^l]Jl2]| , because from our point of view, 
A tg is composed mainly of the temperature independent 
Apg. At optimal doping the above ARPES experiments 
show a significant change of the excitation gap itself with 
decreasing temperature in addition to the sharpening of 
the peak. From our point of view, this significant change 
of A t g (at optimal doping) with temperature is due to 



To find out an empirical relation between A cg , A tg and 
Apg, we plotted i/A| g — Apg (long-dashed curve) and 



A cg near the transition. 



Atg — Apg (dashed curve). The former is nearly indepen- 
dent of doping, in agreement with the result often found 

by other authors using A C5 = yj &ig ~ ^ P g @ to extract 
the coherence gap from experiment. Within our numeri- 
cal uncertainties, the coherence-gap A cg is approximately 
given as A tg — A pg and its behavior is easier to recon- 
cile with the doping dependence of T c and d. Recently 
Krasnov H] also found that the coherence gap (dubbed 
as the SC gap in his paper) shows a doping dependence 
similar to that of T c . 



Since the zero-temperature critical field H C 2 is a prop- 
erty of the condensate, the above V e ff model can also 
be used to estimate it. Using the d-wave estimate pq] 



(ttA, 



and 



H c2 = 0.521$ (2tt£o) with £ = Hv F ; 
$o = hc(2e)' 1 , we find that H c2 (x = 0.08) = 42T and 
H C 2 (% = 0.15) = 151T. The experimental values ]5{l are, 
respectively, around 30T and 60T, somewhat smaller. 
However, there is still disagreement on the experimen- 
tal value of H C 2 and recent estimates are larger |R| . The 
above model is extremely crude since it does not take 
into account the anisotropy and interaction renormaliza- 
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tion of the Fermi velocity. Within this crude model, we 
have that H C 2 will decrease at dopings larger than opti- 
mal because above optimal doping A cg decreases. 

VI. DISCUSSION 

In this Section we compare the present results with the 
predictions of some of the leading theories for the high 
T c superconductors and we provide further comparison 
with experiments. Before we proceed, we should stress 
that the starting point of our approach is meant to apply 
to the intermediate to strong-coupling regime and thus 
should be distinguished from spin-fluctuations theories 
that apply from weak to intermediate coupling. fil| . |62f 

First, consider the question of why the superconduct- 
ing T c decreases near half-filling. As is clear by now, 
in our study the reduction of T c near half-filling occurs 
already at mean-field level due to the competition with 
AF correlations (or local moment formation) that open 
up a mean-field gap, making SC order less favorable. Al- 
though SC fluctuations will indeed make T c decrease with 
decreasing doping because of the smaller superfluid den- 
sity |^,^], such fluctuations are not necessary to have 
the correct qualitative behavior. 

Second, what is the nature of the normal state pseudo- 
gap? In this paper the pseudogap is just a crossover phe- 
nomenon due to short-range equal-time spin correlations 
which induce short-range dynamical fluctuations in the 
particle-particle channel. The latter are very effective in 
creating a pseudogap for several reasons. First note that 
the mass renormalization makes the bandwidth and the 
Fermi velocity quite small. In addition, the pseudogap 
opens up near the band edges where the local Fermi ve- 
locity is smallest. This means that the correlation length 
associated with the fluctuations does not need to be large 
to be in the renormalized classical regime. Second, for 
q = d-wave singlet fluctuations, no umklapp condition 
needs to be satisfied contrary to antiferromagnetic hot 
spots, f&ifl The proposed origin for the pseudogap con- 
trasts with some of the recent proposals in which mean- 
field flux phase, p5| circulating current phase, [ |r36| and d- 
density wave phase |p6|j are interpreted as the pseudogap 
state. In these scenarios, the pseudogap state has bro- 
ken symmetry such as time reversal, translational and 
rotational symmetries, and the crossover-like behavior 
observed in experiments is argued to be caused by an 
impurity effect. 

Third, let us consider experimental evidence that there 
are significant short-range AF correlations below x c . Re- 
cent muon spin relaxation 67 6^] and ac-susceptibility 
measurements, show the existence of slow spin fluctua- 
tions below x = 0.19. By using polarized and unpolarized 
elastic neutron scattering experiments as well as zero- 
field muon spin resonance, Sidis et al. p9] | observed an 
unusual commensurate AF phase on a nanosecond time 



TABLE I. Temperature dependent antiferromagnetic cor- 
relation length evaluated for two dopings using fluctuation 
theory. The fluctuations may either be commensurate (C) or 
incommensurate (IC). The mesh used in the calculation was 
128 x 128 hence the largest number in this table is a rough 
estimate. 





x = 0.1 


x = 0.15 


T = 0.02J 


133 (IC) 


4.1 (IC) 


T = 0.04J 


25.5 (C) 


2.6 (IC) 


T = 0.06J 


6.1 (C) 


1.9 (IC) 



scale that coexists with superconductivity in underdoped 
YBa 2 Cu306.5- Lake et al. |70| found AF correlations in- 
side a vortex core in underdoped cuprates. Most recently, 
Hodges et al. [fnj have even shown that the addition of Co 
to an optimally doped YBCO compound induces anti- 
ferromagnetism above T c that survives and coexists with 
d — wave superconductivity below T c . The coexistence 
seems to be at the microscopic level with £ > 200A. In 
the present paper, the presence of AF correlations be- 
low x c does not come from any extra symmetry between 
AF and SC correlations such as SO(5) symmetry, but 
from the fact that the J part of the t — J Hamiltonian 
has those competing correlations in the low energy sector. 
Table 1 gives a few results for the temperature-dependent 
antiferromagnetic correlation length in our approach for 
two dopings. As expected, the correlations decrease with 
doping. An upper bound, order of magnitude estimate 
for the value of the static local moment is given by the 
value of m in Fig.|l](b) 

Fourth, we point out that there is numerical and exper- 
imental evidence that short-range spin correlations are 
closely related to the normal state pseudogap. In their 
recent calculations obtained from the dynamical cluster 
approximation (DCA) for the Hubbard model, Jarrell et 
al. H noted that the fall of T* with doping is closely 
tied to the diminishing of the short-range spin correla- 
tions. In the exact diagonalization study of the t — J 
model, Sakai and Takahashi |72| found that the pseudo- 
gap behavior is associated with the development of static 
AF spin correlations with decreasing temperature. In 
his high temperature series expansion study of the t — J 
model, Putikka (73) showed that pseudogap crossover oc- 
curs when static AF correlation length is about one lat- 
tice spacing. By analyzing various experiments (Raman, 
spin- lattice relaxation rate, ARPES, Zn-substitution ef- 
fects, inelastic neutron scattering experiment and so on), 
Tallon and Loram |l8| concluded that the pseudogap is 
intimately connected with short-range AF correlations. 
Short-range dynamical AF correlations can also cause 
the normal state pseudogap through AF spin fluctua- 
tions. As explained above however, it was found in a 
recent paper |4(| that the pseudogap always appears ear- 
lier in the dynamical particle-particle channel than in the 
dynamical AF spin fluctuation channel, even though the 
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short-range equal-time correlations that generate these 
dynamical fluctuations are due to AF correlations. One 
of the difficulties with AF spin fluctuation scenario of the 
pseudogap is how the normal state pseudogap happens 
to have the same momentum dependence and magnitude 
as the SC gap near T c , J74| as ARPES and tunneling ex- 
periments show. In our case, this occurs naturally. Our 
results remain qualitatively the same when t' is —0.34. 

Fifth, is x c ~ 0.19 — 0.20 accidental for the special 
choice of t/J = 3? According to the calculations of the 
slave-boson mean-field solution of the Hubbard model 
|75|[76| , the onset of short range AF correlations starts 
to appear from x = 0.20 — 0.21 for a wide range of U = 
7t — 16t (See Fig.l of Ref. J76| for more details). In our 
approach, x c varies from 0.22 to 0.17 when t/J varies 
from 2.5 to 4. It is believed that a realistic strength of 
the Coulomb repulsion U = 4t/ J is of the order of the 
bare bandwidth 8t in two dimension. j77j 

Last, through the microscopic separation of hole-rich 
SC regions from AF regions, ]7§] the stripe structure 
tends to maintain AF correlations more effectively than 
the other case in which they are uniformly suppressed 
by doped holes. Then one can surmise that in the stripe 
state T c is somewhat suppressed from that in the uniform 
state. While the differences between the stripe state and 
the uniform state can be quite substantial, we do not 
see the stripes as necessary to obtain superconductivity. 
The main results reached in this paper are not expected 
to qualitatively change in the presence of a dynamical 
stripe structure. 

VII. CONCLUSION 

In this paper we proposed a simple phenomenological 
procedure that allows one to study the competition be- 
tween antiferromagnetism and <i-wave superconductivity 
in the high T c cuprates. The no-double-occupancy con- 
straint is taken into account in the effective mass approx- 
imation. Correlation functions that are local in space and 
time are evaluated from a mean-field factorization of the 
t — J Hamiltonian in both the AF and the d— wave SC 
channels. (The d— density wave channel does not con- 
tribute in the whole doping range.) The local correlation 
functions are then used in the fluctuation-dissipation the- 
orem to compute renormalized vertices that allow one 
to obtain the full dynamical susceptibilities. In other 
words, in this approach, equal-time correlation functions 
determine the value of effective vertices in all available 
channels. In particular, short-range equal-time AF cor- 
relations determine the value of effective interactions 
in the particle-particle channel. This effective weak to 
intermediate-coupling approach cannot work for x < 0.08 
where the bandwidth becomes less than the interaction 
strength. 



The calculated total excitation gap A tg (in the su- 
perconducting state) and the normal state pseudogap 
A pg are in good agreement with experimental results, 
as shown earlier in Ref. E(J . Obtaining a description of 
superconducting properties arising from a highly corre- 
lated state remains a challenge. Here we have taken the 
simplest approach. We compute an effective strength 
of pairing interactions which takes into account (pair 
breaking) AF correlations, V e ff, by requiring that the 
zero-temperature order parameter d obtained with V e ff 
equals that obtained in the full mean field equations cor- 
responding to Eq.(^). This effective interaction allows 
one to obtain properties of the condensate (and of the 
condensate only), namely the condensation energy U(0), 
the coherence gap A cg that has been observed in Andreev 
reflection (56) as well as H c i- The calculated coherence 
gap closely follows the doping dependence of T c or d, 
and is approximately given as A cff ~ A tg — A pg within 
our numerical uncertainties. A qualitative phase diagram 
for the cuprates may thus be obtained. The systematic 
decrease of U(0), A cg and H c2 with decreasing doping 
below x c ~ 0.2 can be understood as a result of the com- 
petition between AF and SC order occurring in the low 
energy sector of the mean-field t — J Hamiltonian. On 
the other hand, in the present picture the systematic in- 
crease of Ap g and A tg when x decreases below x c ~ 0.2 
is due to the growing short-range equal-time spin corre- 
lations which induce dynamical singlet fluctuations with 
(i-wave symmetry, low characteristic frequency but with 
small correlation length. As pointed out before in Ref. 
p0| one prediction of this approach is that the pseudogap 
should survive if the superconducting state is destroyed 
by a magnetic field. This is consistent with the obser- 
vation of antiferromagnetic ordering inside vortex cores. 
j70f Neutron scattering experiments should be able to 
check from energy integrated structure factor that the 
short-range spin correlation functions are in good agree- 
ment with the mean-field predictions. Finally, we expect 
H C 2 to decrease with doping in the overdoped region. 

Although some of the results of this paper are in quan- 
titative agreement with experiment, the overall descrip- 
tion is definitely qualitative only. The main disagreement 
with experiment is that the renormalized Fermi velocity 
that we find is three to four times smaller than measured. 
A satisfactory theory of high T c will necessitate a better 
treatment of strong-coupling short-range correlations. 
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